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Consider a thermodynamic system composed of two independent subsystems, A and 
B, in contact to each other. The concept of "composability" [1] states that total entropy 
of any kind, S (A, B), is given by a function of the same type of entropies of the 
subsystems, S(A) and S(B): 



S{A,B) = f{S(A),S(B)). (1) 



Here,/is a certain bivariate function of the C 2 -class and is assumed to be symmetric 



f{S(A),S(B)) = f{S(B),S(A)). (2) 



The celebrated Boltzmann-Shannon entropy possesses the additivity property [2] 



f{S(A), S(B)) = S(A) + S(B). (3) 



Another example may be supplied by the rule 



f(S q (A), S q (B)) = S q (A) + S q (B) + Q(q) S q (A)S q (B), (4) 



where Q(q) is a function of the "entropic eindex", q, satisfying 2(1) = 0. 
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Q(q) = l-q and Q(q) = q-l correspond to the Tsallis entropy [1] and the modified 
Tsallis entropy [3], respectively. In both cases, the derivation of q from unity measures 
the degree of nonextensivity of the entropies. Clearly, additivity in Eq. (3) is recovered 
in the limit g— >1. The property in Eq. (4) is referred to here as Tsallis-type 
pseudoadditivity. 

To construct thermodynamics based on entropy, it is necessary to define the state of 
equilibrium, first. This is essentially relevant to the zeroth law of thermodynamics. In 
this Letter, we derive the most general form of pseudoadditivity of composable entropy 
based only on existence of equilibrium. 

Consider certain extensive variables, X(A) and X(B), of the subsystems, A and B, 
respectively. X may be, for example, internal energy, system volume or particle 
number. Then, the equilibrium state [4] characterized by the maximum of the total 
entropy in Eq. (1) with the fixed total amount, X(A, B) = X(A) + X(B) , yields 

df{S(A),S(B)) dS(A) = df{S(A),S(B)) dS(B) 
dS(A) dX(A) dS(B) dX(B)' 

To establish the zeroth law of thermodynamics, by which physical variable (such as 
temperature, pressure or chemical potential) common to the subsystems in equilibrium 
can be identified, it is essential to be able to realize factorization of Eq. (5) in the form 
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F(A) = F(B). 



(6) 



This means that the following set of equations should hold: 



df(S(A), S(B))_ = k{s{Al S{B)) g{s{A)) (7) 



dS(A) 

df{S(A),S(B)) 
dS(B) 



= k{S(A),S(B))h{S(A))g(S(B)). (8) 



Here, g and h are some functions, and k is a differentiable bivariate function. In 
particular, h has to be differentiable. The function, k, does not have the factorized form, 
in general. The symmetry off shown in Eq. (2) tells us that k is also symmetric 



k(S(A),S(B)) = k(S(B),S(A)). (9) 



Since/is of the C 2 -class, the integrability condition holds, leading to 



dJ ^Wr l hm) S{S(B)) +k{s(Al m dJ Si S{S{B)) 

Using Eqs. (7) and (8), we rewrite this equation as follows: 



1 dk{S(A),S(B))df{S(A),S(B)) [ dh{S(A)) 



k 2 (S(A),S(B)) dS(A) dS(B) dS(A) 



1 dk{S(A),S(B))df{S(A),S(B)) | s ^ A ^ dh{S(B)) (n) 



k 2 (S(A), S(B)) dS(B) dS{A) v " dS(B) 

This is an identity and therefore leads to the following equations for the factorized and 
non-factorized parts: 



dk{S(A), S(B)) df{S(A), S(B)) _ dk{S(A), S(B)) df{S(A), S(B)) 

dS(A) dS(B) dS(B) dS(A) ' 1 



respectively. The general solution of Eq. (13) is given by 



k{S(A), S(B)) = G(f{S(A), S(B))), (14) 



where G is an arbitrary differentiable function. 

First, we consider the simplest case when A: is a constant function. Without loss of 
generality, such a constant can be set equal to unity. Then, Eqs. (7) and (8) become 
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"^^--^M (16) 



respectively. From Eq. (12), it follows that 



1 dh(S(A)) _ 1 dh{S(B)) =x 



g{S(A)) dS(A) g{S(B)) dS(B) 



where X is a separation constant. Using Eq. (17), we rewrite Eqs. (15) and (16) as 



*« = i«i ilM , (18) 



Here, the case A = has to be interpreted as the limit X — > . Integrating these 
equations, we find 



f x (S(A), S(B)) = jh x {S(A)) h x {S(B)) + const . (20) 



To have the convergent result in the limit X — > , without loss of generality, we can set 
the constant term in Eq. (20) equal to -Ilk and simultaneously impose the condition 



\imhAS) = \, (21) 



for V S. Thus, we obtain 



/lWW= ««h! , (22) 



This is the form of the function of composability prescribed by existence of equilibrium. 

If both of the subsystems are at the completely ordered states, then S(A, B) = 0. 
This means that 



M0) = 1, (23) 



for VA . On the other hand, if only the subsystem B is at the completely ordered state, 
then we have 



S{A,B) = S(A). (24) 
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Therefore, from Eq. (22), we conclude 



5(A), 



(25) 



or equivalently, 



h x (S) = l + XS. 



(26) 



With this form, Eq. (21) and (23) are clearly fulfilled. Thus, we find that Eq. (22) is 
equivalent to Eq. (4) with the identification, X-Q. In other words, Tsallis-type 
pseudoadditivity corresponds to the simplest case when k in Eqs. (7) and (8) is a 
constant function. 

In connection with the above result, we wish to mention that in Ref. [5] a 
nonextensive generalization of the Shannon- Khinchin set of axioms for the Boltzmann- 
Shannon entropy is given. The axioms presented there are 

[I] S (p v p 2 ,L , p w ) is continuous with respect to all its arguments and 



takes its maximum for the equiprobability distribution p ; = 1 / W 



(i = l,2,L ,W), 



[II] S q (A,B) = S q (A) + S q (B\A) + (l-q)S q (A)S q (B\A), 
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[Ill] S q (p l3 p 2 , L , p w , p w+l = 0) = S q (p x , p 2 , L ,p w ). 
Here, p ; (z = 1, 2, L , is the probability of finding the total system in its zth state. 
S q [B\A) is the conditional nonextensive entropy [5] of the subsystem, B, given the 
sunsystem , A. Comparing this set with the Shannon-Khinchin one, the only difference 
between the two is in [II]. (The Shannon-Khinchin axioms are recovered from [I] -[III] 
in the limit q—>l.) In a particular case when A and B are independent each other, 
S q [B\A) = S q (B) holds, and accordingly [II] becomes the Tsallis-type pseudoadditivity 
rule in Eq. (4) with Q(q) = l-q. The uniqueness theorem proved in Ref. [5] states that 
a quantity S q satisfying [I] -[III] is equal to the Tsallis entropy [1] 



s q {Pi'P2> L >Pw) = r— 
1 - q 



EW- 1 



(27) 



with q>0. 

Next, let us discuss the general case in Eq. (14) and put 



a W-d5W <28) 

df 



where H is a certain differentiable function. Then, Eq. (22) is now replaced by 



9 



^ S) )=«»)H. (29) 



Here, we are using the same notation as in Eq. (22), which may not cause any 
confusions. For the completely ordered subsystems, we have 



«(0)= A %i. ( 30) 



Also, if only the subsystem B is at the completely ordered state, then we have 



H(S(A))= h ^\ h ^-\ (31) 
A 



or equivalently, 



MSM) )_!±*M 

MO) 



Therefore, from Eqs. (29) and (32), we find 



H(s(A B)] = H{S(A)) + //(S(*)) + A H(S(A)) H(S(B)) - H(0) 
y y ' >> 1 + Aff(0) 
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In a particular case when H(0) can be taken to be zero, this equation becomes 



H(S{A, B)) = H(S(A)) + H(S(B)) + X H(S(A)) H(S (B)). (34) 

Eq. (33) is the most general pseudoadditivity rule for entropy prescribed by existence of 
equilibrium. Tsallis-type pseudoadditivity is recovered when H is the identity function. 

As an example, let us consider the choice: H(f) = and X — > . In this case, Eq. 
(34) gives S(A, B) = S(A) + S(B) + 2^S(A)S(B) . This rule might be relevant to the 
black hole entropy, which is proportional to the square of the mass of the black hole [6]. 

In conclusion, we have derived the most general pseudoadditivity rule for composite 
entropy based only on existence of equilibrium. We have shown how Tsallis-type 
pseudoadditivity can be obtained as the simplest case. 

This work was supported by the Grant-in-Aid for Scientific Research of Japan 
Society for the Promotion of Science. 
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